Persistent currents in mesoscopic cavities with surface roughness 
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0^ . We consider a bidimensional discrete annular cavity with surface roughness (SR) threaded by 

a magnetic flux. In the ballistic regime and at half filling, localized border-states show up. These 
^ \ border states contribute coherentely to the persistent current and the magnitude of the typical 

O ■ current Ityp is enhaced with respect to their value in the absence of confinement. This result is 

^^ ' robust with respect to changes in the relevant lengths of the sample. 



1-^ . In recent years, advances in nanotechnology made possible to design mesoscopic samples in which the carriers are 

C/j ■ mainly scattered by the boundaries of the system |I|. In this situation, where the elastic mean free path is much larger 

, than the system size, the sample can be considered in the ballistic regime in order to compute the relevant transport 
H ■ quantities. In addition, the potential used to confine the carriers into the mesoscopic device has a symmetry quite 

1 t [ different from that of the lattice of the sample. Therefore, even when the bulk disorder is absent, surface roughness 
Ci i (SR) is present with the consequence of the irregular scattering of the carriers at the borders of the sample. 

H ' In this article we present a model of SR which emerges naturally when the topology of the discrete lattice and the 

I . confining potential are different. The relevant feature of our model is that the effect of the SR on the eigenstates 

depends on the energy of these states. As the energy spectrum is bounded, the lattice parameter a is a cut off for the 

allowed wavelengths and the wavelenghts close to a will be the most sensitives to the SR. Therefore, the SR that we 

fj ' model will be relevant near half filling, i. e in the metallic regime. 

1—^ [ When the sample is threaded by a magnetic flux persistent currents are originated P,p| . Although much work has 

, ' been done concerning the influence of bulk disorder on persistent currents 0, less is known about the ballistic regime. 

►^ Our goal is to characterize the effect of the confinement and the SR on the magnitude of the persistent current in 

r — i that regime. 

p^ The total current can be calculated at zero temperature as / ^ dE/d(j>, with E the total energy of the system and 
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(f> the magnetic flux. The value of Ityp — \/ J P d(f> depends on the number of open channels present in the system 

M. For discrete sytems with cylindrical geometry and away from half filling, Ityp ~ VM ||,^. This behavior can be 
U^ understood taking into account the properties of the crossings (or quasicrossings) that appear in the spectrum when 
^^ ■ the flux is varied M. They occur mainly between levels belonging to different channels and frustrate the coherent 
increased of the current B . Although very particular geometric configurations can give rise to coherent contributions 
of all the channels MM , those sharp configurations are very difficult to obtain experimentally. 

We consider a. N x N cluster on a square lattice. In the following all the lengths will be given in units of a. Taking 
i-pj . a single atomic level per lattice site the Hamiltonian is, 

O. H ^^em,n\m,n><in,n\+ ^ tjn,n;m',n'\m,n >< m' ,n'\ , (1) 



m,n:m' ,n' 



• ^ . where (to, n) = m labels the coordinates of the sites in the lattice. The hopping integrals tm;m' are restricted 

^\ ' to nearest neighbours. Assuming that the vector potential A has only azimutal component we take tm;m+i = 

C^ ' t exp (i J A.dl) , in units of the quantum flux (f>o, where 1 is a vector that points from the site m to any of its 

four nearest neighbours. 

With the on-site energies, em, we simulate the conflning potential in order to obtain the required profile of the sample. 

In particular for an annular cavity with internal radius r and external radius R, we define 

ifr<|TO|<i? ,. 

V otherwise , ^ ' 

where V is large enough to confine the carriers within the annulus. For the present geometry the two relevant lenght 
scales can be chosen as the internal radius r and the mean width of the sample W. While the second determines 



M for a given value of the Fermi energy jl^l, the ratio s = W/r controls the strength of SR. For s >> 1 and for 
carriers with a wavelength of the order of a, two succesive scattering events by the inner and outer boundaries of 
the sample are uncorrelated. On the other hand, for s << 1 and the same wavelengths those events are strongly 
correlated. Therefore the eigenstates corresponding to these particular wavelenghts will evolve from generic extended 
ones to strongly localized as the value of s is reduced. 

In order to describe our results we will focus in two representative samples: SI, with W — AQ and r = 20 and S2, 
with W — 16 and r — 90. The values are given in units of a. In both samples the averaged area of the sample is A 
and equal to the total number of sites. 

In Fig. Q^ we show for sample 5*1 the bottom (Fig. 0(a)) and the quarter filling (Fig. mh)) regions of the energy 
spectrum as a function of the normalized flux a = (t>/4'o together with a generic charge distribution for each region. 
Fig. p](a) corresponds to the large wavelength numbers, that are the less sensitive ones to the effect of the SR. This 
region of the spectrum looks qualitative similar to the spectrum of an integrable Aharonov-Bohm annular billiard [p| . 
That is, crossings between levels belonging to the same channel at a = and a = 1/2 and crossings between levels 
belonging to different channels for < a < 1/2. Moreover, a generic charge distribution corresponds to eigenstates 
with well defined value of the angular quantum number. 

In the quarter filling region of the spectrum the typical wavelengths are greater than a but small enough to make 
the states sensitive to the global polygonal shape of the sample. Therefore, as the energy increases, avoided crossings 
are showed up as a function of a. The corresponding charge distributions are quite similar to those of a continuous 
polygonal billiard and are extended states in general (see Fig. 1(b)). In both energy regions, when an occuppied 
level crosses (or quasicrosses) an empty one, the total current / exhibits a discontinuity (or abrupt oscillation). The 
magnitude of Ityp will be reduced with respect to the situation with no crossings in the interval (0,a/2) y,|[. As 
it is shown in Fig. H, below half filling, Ityp is a highly fluctuating function of the number of particles Np. This is 
a consequence of the randomness in the distribution of crossings (or quasicrossings) present in these regions of the 
spectrum 18]. For sample 52 a qualitative similar behavior is found. 

The upper region of the spectrum, close to half filling, is shown in Fig. g(a) for 5*1. As we have already mentioned, 
in this region the states are sensitive to the fine structure of the boundaries of the sample. A bunching of quasi- 
degenerated eigenenergies appears at zero fiux corresponding to two different types of localized states. One type are 
the angular localized states whose characteristic charge distribution is depicted in the left side of Fig. ||(b). These 
states are trapped in regions where the sample has more steps on the boundaries. They are unsensitive to the magnetic 
flux and are the origin of the flat lines in the spectrum. The other type of states are mostly radially localized and 
look like whispering gallery modes. These states are sensitive to the magnetic flux and therefore will carry a finite 
current. The right side of Fig. ||(b) displays the characteristic charge distributions for these states. We remark 
that qualitatively the same description holds for sample S2. Nevertheless, as the fraction of angular localized states 
decreases with the width W of the sample, the number of quasidegenerated fiat levels in the spectrum of S2 at half 
filling will be greater than in 5*1. 

Opposite to what happens for models with diagonal disorder, for the present system the localization is obtained 
without introducing an additional scale of energy. Only remains the scale defined by the kinetic energy. Therefore, 
at zero fiux, the energy of the localized states will be very close to zero. This is a consequence of the finite character 
of the spectrum and its symmetry. 

In order to quantify the described effect we have to consider two energy scales. One is the mean level spacing A 
near half filling between levels of a given channel and in the absence of confinement, A = — , where p is the mean 
density of states near half filling B. The other relevant energy scale, A^, is a measure of the effect of the confinement 
on the eigenstates and eigenvalues of the system without borders. To obtain a generic quantitative relation between 
Ag and A is not a simple task. However we present the following argument: for W —> oo the effect of the confinement 
tends to zero, and then Ag < A/M. On the other hand, ior W ^ the effect of the SR is maximized and Ag > A. 
In this regime the eigenstates become angular localized, i.e are not current carrying. By increasing W we can go 
from one limit to the other one. Therefore in a wide intermediate regime it should be satisfied that ^ < Ag < A. 
When this condition is acomplished the interaction with the borders mixes states belonging almost all to different 
channels (i.e mainly with similar angular momentum), making radially localized eigenfunctions, which corresponds 
to the quasidegenerated levels bunching at zero energy that are current carrying. Therefore even in the presence of 
SR, typical states like the whispering gallery modes survive in the discrete system close to half filling. This have an 
important consequence on the values of Ityp close to half filling, as it is shown in Fig. for both samples SI and S2. 
As there are not crossings between occupied and empty levels in all the region of variation of the flux, Ityp increases 
monotonically with Np. A plateau occurs at the maximum value of hyp, and appears when the flat levels (angular 
localized states) begin to be filled. As we have already mentioned the width of the plateau for 51 (Fig. 0(a)) is smaller 
than for 52 (Fig. |(b)). 

As far as we know the enhacement of the persistent currents due to confinement effects in discrete systems has 



not been adressed in previous theoretical works. We believe that the present work could help to understand early 
experimental resuls in Au rings [[l3[ which have shown big discrepancies with the existing theories. 
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FIG. 1. (a) Energy levels at the bottom of the spectrum as a function of the rescaled magnetic flux a for sample SI together 
with a typical charge distribution, (b) Idem as (a) but for the quarter filling region of the spectrum. 
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FIG. 2. (a) Energy levels close to half filling as a function of the rescaled magnetic flux a for S'l.(b) Typical charge distri- 
butions for angular localized states (left side) and radially localized states (right side). 
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FIG. 3. (a) Typical current as a function of A'p for sample 5*1. (b) Idem as (a) for sample S2. 



